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Abstract. We review (i) observations and numerical simulations of vortical flows 
in the solar atmosphere and (ii) measurements of the horizontal magnetic field in quiet 
Sun regions. First, we discuss various manifestations of vortical flows and emphasize 
the role of magnetic fields in mediating swirling motion created near the solar surface 
to the higher layers of the photosphere and to the chromosphere. We reexamine exist- 
ing simulation runs of solar surface magnetoconvection with regard to vortical flows 
^ ^ and compare to previously obtained results. Second, we review contradictory results 

and problems associated with measuring the angular distribution of the magnetic field 
in quiet Sun regions. Furthermore, we review the Stokes-V-amplitude ratio method for 
the lines Fe i AA 630.15 and 630.25 nm. We come to the conclusion that the recently 
discovered two distinct populations in scatter plots of this ratio must not bee interpreted 
in terms of "uncollapsed" and "collapsed" fields but stem from weak granular mag- 
netic fields and weak canopy fields located at the boundaries between granules and the 
intergranular space. Based on new simulation runs, we reaffirm earlier findings of a 
predominance of the horizontal field components over the vertical one, particularly in 
the upper photosphere and at the base of the chromosphere. 



1. Introduction 

Two different topics are reviewed in this contribution. (1) Vortical flows, which have 
been observed to exist in the deep photosphere and in the chromosphere. Over the 
past three years, they have anew attracted the attention from observers and computa- 
tional physicists. We attempt to find and emphasize interrelations between the various 
manifestations of vortical flows in the solar atmosphere. (2) The polarimetry of the hor- 
izontal magnetic field in quiet Sun regions. We review the disparate results that have 
been obtained in the past and highlight the difficulties associated with such measure- 
ments. We discuss biases and propose an alternative method for dealing with selection 
criteria. The main part of the second chapter is devoted to the Stokes-V amplitude-ratio 
method and the origin of the dichotomy observed in scatter plots of the V amplitude of 
Fe I 630. 15 nm vs. the V amplitude of Fe i 630.25 nm. 



2. Vortical Flows and Vortex Tubes 



Vortical flows, swirls, whirlpools, vortex tubes in the solar atmosphere have become a 
focus of intense research in the past three years. This chapter presents a brief review of 
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this research focus and attempts a comparison of som e of these results with simulations 
carried out with the C05B0LD (Freytag et al."2012') code by the authors. 

Vortical flows in the photosphere of the Sun were reported long before the present 
revival. Most notably is the vortex flow of iBrandt et al.l (Il988h . which was a vortical 
movement of granules that persisted for 1.5 h. These authors conjectured that if such 
vortices were a common feature of the solar convection zone, they might "provide an 
important mechanism for heating of stellar chromospheres and coronae by twisting the 
footpoints of magnetic flux tubes''.^ 

From numerical simulations, iNordlundl (1 19851 ) reports "a circular motion around 
the center of the downdraft, and the circular velocity is amplified as the downdraft 
narrows ('bath-tub' or 'inverted tornado')". The centripetal force associated with this 
vortical flow was so strong that the gas pressure gradient opposing it, apparently led 
to numerical difficulties because the ga s pressure in the ce nter of the downdraft be- 
came too small or negative. Similar to IBrandt et al. (119881) . he also wondered about 
the consequences of such 'inverted tornadoes' for the higher layers of the atmosphere 
and conjectured that "the fact that any vertical magnetic field lines in the surrounding 
photosphere must be carried towards, and 'sucked into', these downdrafts also makes 
the phenomenon potentially very important as source of hydromagnetic disturbances". 

IWang et all 31995) computed from a time series of continuum images from the 
Pic du Midi Observatory the horizontal vector flow field and associated divergence and 
the vertical component of the curl and found that the curl is correlated with regions of 
negative diver gence, suggesting exce ss vorticity in intergranular lanes. A similar result 
is reported by iPotzi & Brandt! (l2007h . 



2.1. Vortical Flows in the Photosphere 

With the ever increasing spatial resolution and homogeneity of long duration image 
sequences thanks to larger apertures of ground based solar telescopes, to high order 
adaptive optics systems and advanced image restoration techniques, and thanks to space 
based and balloon borne telescopes, vortical flows on ever smaller scales became visi- 
ble. 



Table 1 . Vortex properties as detected bv' Bonet et al.l (120081 120101) . The numbers 
in parentheses in the second column refer to the number of vortices with clockwise 
(first) and counterclockwise (second) sense of rotation. 



telescope 


number of detected 


space-time density 


mean life-time 




vortices [-] 


[Mm^^ minute"'] 


[minutes] 


SST 


138 (68,70) 


1.8 X 10-^ 


5.1 


Sunrise 


42 (15,27) 


3.1 X 10-^ 


7.9 



Thus, iBonet et al.l (l2008h discovered small whirlpools on the solar surface with a 
size similar to terrestrial hurricanes (< 0.5 Mm), using the Swedish Solar Telescope 
(SST). They detected them because some magnetic bright points follow a logarithmic 
spiral on their way to being engulfed by a downdraft. Their image sequences at disk- 
center show 0.9 X 10~^ vortices per square Mm, with a lifetime of the order of 5 minutes, 
and with no preferred sense of rotation (see Table[T]). They repeated this kind of analysis 
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with data from the ballon borne solar telescope Sunrise (lBartholll201ll) and obtained 
similar results (see Table [T]) with the notable difference that from the Sunrise data, 
they detected a preferred sense of rotation of the vortices. They speculate that this 
might be due to the fact that the Sunrise image sequences were recorded away from the 
solar equator, different from the previously obtained SST data. Away from the equator, 
differential rotation may have introduced a preferred sense of rotation. It remains to be 
expl ored if this effect is indeed takin g place. 

IVargas D omfnguez et al.l ( 201lh st udy a region of 69" x 69" of quiet Sun that 
includes the field of view investigated by iBonet et all (12008.) and derive statistical prop- 
erties of s wirl motions in the photosphere. They confirm the value for the space-time 
density of Bon et et all ( 2008h and find no significant preference in the sense of rotation. 

On a slightly larger scale, lAttie et all (2009,) identify two long lasting vortex flows 
located at supergranular junctions, based on two time series of granulation from the 
Hinode/SOT broad band imager. The first vortex flow lasts at least 1 h and is « 20" 
wide, the second vortex flow lasts more than 2 h and is ~ 21" wide. In one case, 
the corresponding magnetogram shows a magnetic element of polarity opposite to the 
magnetic field in the center of the vortex, moving around the vortex center. At the same 
time, there is enhanced Ca u emission at the vortex center, which suggests possible 
magnetic diss i pati on because of entwinement of the magnetic fields. 

iBalmaceda et al.l (1201 Oi) report an observation of magnetic flux concentrations be- 
ing dragged towards the center of a convective vortex motion in the solar photosphere 
from data, including magnetogram sequences, simultaneously acquired with the SST 
and the Solar Optical Telescope (SOT) of Hinode. They argue that these small-scale 
motions were "likely to play a role in heating the up per solar atmosphere by twist- 
ing magnetic flux tubes." Likewise, iManso Sainz et al. (2011 ) find two instances of 
magnetic elements which describe a vortical motion on a scale of < 400 km, using 
spectro-polarimetric and magnetographic data obtained with Hinode. 



2.2. Swirls in the Chromosphere 



Using the Swedish Solar Telescope (SST) with the CRISP Fabry Perot system. IWedemeyer-Bohm & Rouppe van der 



(120091) found in narrow -band filtergrams with a FWHF of 11.1 pm in the line core of 
Ca I 854.2 nm swirling motions in chromospheric layers of a quiet Sun region inside 
a coronal hole. There, the view is not obstructed by chromospheric filaments, which 
gives free sight to chromospheric layers that would not be easily accessible in more 
active regions (see Figure [T|). 

The size of the swirl events is comparable to large terrestrial typhoons but their 
physical origin is a very different one. This is obvious as co-temporal broad band im- 
ages show close groups of photospheric bright points that move with respect to each 
other, right beneath the swirls. Since such bright points are a reliable signature of 
tiny magnetic flux concentrations in the photosphere, it is highly probable that the 
origin of the chromospheric swirls are to be found in these moving magnetic flux 
concentrations. The swirls also exhibit Doppler shifts of -2 to -4 kms"^ with peak 
values of up to -7 kms~^ There is no clear indication of a swirling motion of the 
bright points, but one can imagine that magneti c flux concentratio ns trapped in a pho- 



tospheric swirl similar to the ones detected by iBonet et all (120081 . bOlO ) may start to 



rotate and lead to chromospheric disturbances that would b e similar to the ones ob 



se rved by [ Wedeme ver-Bohm & Rouppe van der VoortI (120091) and similar as predicted 



bv lNordlund(il985l) . 
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Figure 1. Temporal evolution of a swirl event as seen in close-ups of inten- 
sity maps in the wide-band (top row), Ca line wing (upper middle), Ca line core 
(lower middle), and Doppler shift (bottom). The black contours in the bottom 
row mark zero Doppler shift. The grey scale of the Doppler shift is from -5.5 
to -1-5 .5 kms~' with negative values corresponding to blu eshifts and thus upflows. 
Credit: IWedemeyer-Bohm & Rouppe van der Voortl (l2009l) reproduced with permis- 
sion © ESO. 



On a considerably larger scale but also in quiet Sun regions. Izhang & Liul ( 2011 ) 
observe extreme-ultraviolet cyclones in all EUV channels of the Atmospheric Imag- 
ing Assembly (AIA) on board the Solar Dynamics Observatory (SDO). These cyclones 
seem to be rooted in rotating network magnetic fields that are simultaneou sly observed 



in the photosphere with the Helioseismic and Magnetic Imager (HMl; Scho u et al. 



(1201 Ih ). They can last for several to more than 10 hr and are found to be associated with 
EUV brightenings (microflares) and EUV wav es in their later phase. At this stage it is 
unclear if and what kind of connections between lWedemeyer-Bohm & Rouppe van der Voortl ' s 
swirl and iZhang & Lid 's cyclones may exist. 

2.3. Numerical Simulations 



Shelyag et all (120111) carried out magnetic and non-magnetic three-dimensional numer 
ical s imulations of solar granulation with the MURaM code ( Vogleij|201lUVogler et al. 



20051) to analyze the generation of small-scale vortex motions in the solar photosphere. 
Starting with a magnetic field-free model, they observe the generation of swirls in the 
upper part of their model, as soon as a uniform, vertical magnetic field of 200 G strength 
is introduced. They find that the vortices in the upper photosphere are co-spatial with 
the magnetic field concentrations in the intergranular network in the lower photosphere 
that develop a short time after introduction of the magnetic field. Previously already. 
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Vogled (|2004) found such a relation as well. It is of course tempting to identify these 



vortical flows, which obviously oc cur in connection wi th magn etic flux conce ntrations, 
with the chromospheric swirls of IWedemeyer-Bohm & Roup pe van der Voort (2009). 
But for this, maps of Ca i 854.2 nm would need to be synthesized from the simulation 
for a more direct comparison. 

As an independent, qualitative confirmation of the results of IShelyag et al.l ( 2011 ). 
we show in Figure |2] swirling motion that occurs in the upper part of a model at- 



pa 

mosphe re obtained from a simulation with the C05B0LD code (Freytag e t al.l 120121 : 



Beeck et al. 2012) . We have carried out three equivalent simulation runs with a box 



size of 9.6 x 9.6 Mm^ and a depth of 2.8 Mm, where the average depth of = 1 is 
at mid height. Figure |2] (top) shows a time instance of the velocity field as projected 
into the horizontal plane at a height of 1300 km above (t^ ) = 1 for the simulation 
run BO without magnetic fields. Added in gray-scale is the temperature, which shows a 
meshwork of sharp temperature peaks, corresponding to shock fronts. The cellular flow 
field does not correspond to the granular flow field in the deep phot osphere — ^rather it 
is an independent, rapidly evolving pattern at chromospheric heights dWedemeyer et al. 



|2004|) . The longest arrows con^espond to a velocity of about 14 kms"^ No swirls are 
visible in this magnetic field-free simulation. Figure [2] (bottom) shows the correspond- 
ing flow at the same time instance from an equivalent simulation, v50, that started with 
a homogeneous vertical magnetic field of strength 50 G. The flow field is conspicu- 
ously different and the meshwork of shock waves is absent. Swirling, or undulatory 
motion is now omnipresent and must have been introduced through the magnetic field. 
A third equivalent simulation, h50, was carried out in which horizontal field of 50 G 
streng th was advected by updrafts through the bottom boundary, similar to'St einer et al. 
(2008b or lStein & N ordlund (2006). This simulation shows again a flow pattern similar 
to the field-free case shown in Figure [2] (top) because in this case the magnetic field is 
more turbulent and complex and its average strength drops exponentially with height 
because of the many intricate loops that develop. These differences suggest that strong 
swirls in the upper atmosphere may best develop in case of a rather unipolar open field 
configuration as may occur at the base of a coronal hole. 

In another independent study, Moll et al. (2012) carry out comparative simula- 
tions, similar to those shown in Fig. |2] They find a considerably different flow struc- 
ture in the upper photospheric layers of two simulations: the non-magnetic simula- 
tion is dominated by a pattern of moving shock fronts while the magnetic simulation 
shows vertically extended vortices associated with magnetic flux concentrations. They 
note that both kinds of structures induce substantial local heating, in the magnetic case 
thro ugh Ohmic dis s ipatio n associated with the swirling motion. 



Shelyag et al.l (120111) also derive the equation for the vorticity a> = V xvby taking 



the curl of the MHD momentum equation: 



tilting and streching Ti baroclinic term T2 



Dio 1 

— - (C^ • V)V - CJ{V ■ V) + — Vp X V/7gas 

Dt 



magnetic baroclinic term 



magnetic tension 



+ 4vpx 



47r 



+ —V X {B ■ V)B , 

4;rp 



(1) 
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Figure 2. Time instances of the velocity field projected into the horizontal plane 
at 1300 km above (t^) = 1. Top: Magnetic field-free simulation BO. Overplotted in 
gray scale is the temperature. Temperature peaks delineate the meshwork of shock 
fronts. Bottom: Same time instance from the simulation, v50, which started with a 
homogeneous, vertical magnetic field of 50 G. Gray scales in the background show 
the logarithm of the field strength, log |B| from 1 to 100 G. The flow field is strikingly 
different from the field-free simulation and shock waves are absent. Longest arrows 
with a 1% area coverage have v > 14 kms~^ (top) and v > 12 kms~^ (bottom). 
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Figure 3. Vertical component of the vorticity, m~ - (V x v)~ (green and blue for 
opposite signature) together with the surface (brown) of plasma y6 = 1 . At the loca- 
tions of magnetic flux concentrations, the /3 - I surface reaches through funnels into 
the deep photosphere and even below the surface (gray) of optical depth unity. In be- 
tween the funnels, it forms a horizontally extending canopy at chromospheric levels. 
Above the - 1 surface, vorticity is mainly confined to the region of /? << 1 . Below 
it, in the convection zone, vorticity is virtually everywhere present bel ow intergran- 
ular l anes. The simulation was carried ou t with the C0 5B0LD code dFrevtag et al.l 
120 12h . the visualization with VAPOR 3D (IClyne et al.ll2007.) . Courtesy, Ch. Nutto. 



and study the individual sources for vorticity, Ti to T4, in their simulation. Here, v is the 
velocity, p the density, pg^s the gas pressure, and B the magnetic field. They found that 
in the upper photosphere, the magnetic tension term is most important, while beneath 
the surface of optical depth unity, = 1, the baroclinic term \T2\ dominates \T4\. Near 
Tc = ^,\T2\ and \T4\ are of similar magnitude. The other terms are of less importance. 

Since the plasma yS (ratio of thermal to magnetic pressure) generally decreases with 
height and assumes values smaller than one in the chromosphere and within magnetic 
flux concentrations in the photosphere, the results of Shelyag et al. (2011.) indicate that 
vorticity may be generated in regions of small /3 through the magnetic tension term, 
i.e. by Lorentz forces. In fact. Figure |3] supports this conjecture. At the location of 
magnetic flux concentrations, the = I surface (brown) reaches through funnels deep 
into the photosphere and even below the surface (gray) of optical depth unity. Above 
the surface of optical depth unity, vorticity (blue and green) is confined to within these 
funnels where /J « 1. There is little vorticity outside the funnels in the photosphere, 
where yS >> 1. Clearly, this enhanced vorticity within the funnels must be generated by 
the magnetic field, possibly through the vortical motion flux concentrations are forced 
to follow wheii getting t rapped within intergranular swirling sinkholes of the type found 
by iBonet et aL or Wan g et al.1 



Kitiashvili et all (120111) discover in their (field-free) simulations tiny circular 'den- 



sity holes' in horizontal cross sections close to the solar surface. These 'density holes' 
with a density deficit of up to 60% are typically located at vertices of intergranular lanes 
but sometimes they also occur in the middle stretch of an intergraular lane. They are 
associated with a deficit in temperature of about 20%, a downdraft of up to 7 kms~^ 
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and swirling motion with often superso nic horizontal ve locities. iKitiashvili et all s 
"whirlpools" must be the same object that lNordlunj (Il985l) found to be a consequence 
of the centripetal force associated with a vortical fl ow. Underdense n ear-surface ver- 
tical vortices of similar properties are also found bv iMoll et all (1201 ih from magneto- 
convection simulations with the MURaM code ('Vogler 2011 : Vogler et al. '2005'). Ac- 
cording to these authors, they cause a local depression of the optical surface, which 
would potentially be observable as bright po ints in the dark intergranular lanes. 

Most important, IKitiashvili et all (l201ll) also found that whirlpools can attract and 
capture other swirls of opposite vorticity. They showed that this processes of vortex 
interaction, can cause the excitation of acoustic waves on the Sun. They conclude that 
strongly interacting vortices in the top layers of the convection zone, in particular the 
interaction of the vortices with opposite sign of vorticity, may play an important role 
in the excitation of solar acoustic oscillations. In their simulations, they did not find a 
preference in directions of the vortex rotation, which should be expected since neither 
the Coriolis force nor differenti al rotation is in cluded. However, if the findings of a 
preferred sense of rotation of Bo net et al] ( 2010h should substantiate and if theses swirls 
are indeed an important source of acoustic oscillations, one might speculate that this 
could have observable consequences for the p-mode oscillations. One could possibly 
expect reduced power at latitudes of enhanced differential rotation. 

Figure|4](top) shows the bolometric intensity of a simulation snapshot of the above 
mentioned field-free C05B0LD simulation BO. Superimposed on the intensity map is 
the velocity field projected into the horizontal plane at the average optical depth unity. 
This is approximately, what could be derived from observations with the help of correla- 
tion tracking and equivalent techniques. The figure demonstrates, that swirling motions 
are not particularly conspicuous. In most vertices of inter-granular lanes, there is hardly 
any sign of swirling motion visible in this snapshot. Only at (x, y) = (6000, 3000) km, 
there is clear evidence of a 'whirlpool'. The situation changes drastically when ex- 
amining the density at about 100 km below mean optical depth unity, which section 
is shown in Figure |4] (bottom). There, we can see small underdense circular areas 
with a diameter of ~ 100-200 km, which are associated with swirling motion, e.g., at 
{x,y) - (3800, 3200) km, {x,y) = (5000, 8500) km, {x,y) - (1900, 8250) km, and twin 
swirls at (x, y) = (5000, 7000) km. They have a density of typically 60% of the mean 
density, while the intergranular lanes have typically 115%. They also are typically 
25-35% cooler than the average temperature. T hese findings qualitativ ely corroborate 
corresponding results by iKitiashvih et all (|201 ih and lMoU et al.l (l201lh . We do not see 
a bright poi nt in the bo lometric intensity map at the location of the whirlpool as was 
reported by IMoU et all (|2011) but this may be due to insufficient spatial resolution of 
the present simulation. There is no one-to-one correspondence of swirls visible at the 
solar surface and underdense vortices 120 km below it. For example, the conspicuous 
swirl visible at {x, y) = (6000, 3000) km at the surface has no underdense counterpart 
120 km beneath. On the other hand, the underdense swirl at {x,y) = (3800,3200) km 
seems to be associated with an intensity enhancement at the surface but not with a clear 
vortical flow. The twin swirls at (x, y) = (5000, 7000) km are associated with a granular 
lane visible in the intensity map (Fig. |4]top). In fact, it looks like the two 'footpoints' 
of the horseshoe-like granular lane were leading into the twin sinkholes, forming all 
together a loop-like vortex tube. Granular lanes and vortex tubes are reviewed in the 
following subsection. 
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Figure 4. Top: Time instance of the bolometric intensity and the velocity field 
projected into the horizontal plane located at (tc) - 1. Bottom: Corresponding time 
instance of the logarithm of the density (-6.25 < logp < -6.86) and the velocity 
field projected into the horizontal plane located 120 km below (tc) - 1. Under- 
dense vortices ar visible at, e.g., (x,y) = (3800, 3200) km, {x,y) = (5000, 8500) km, 
{x,y) = (1900,8250) km, and twin swirls at (x,y) = (5000,7000) km. They have 
a density of typically 60% of the mean density. Intergranular lanes have typically 
115%. The snapshot is from a magnetic field-free simulation with C05B0LD. 
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2.4. Granular Lanes and Horizontal Vortex Tubes 

It was known for long that in the visible continuum, the edges of granules tend to 
be brighter than the interior of granules or the intergranular space. High resolution 
observations have refined this knowledge. It was found that granules frequently show 
substructure in the form of lanes composed of a leading bright rim and a traiUng dark 
edge, which form at the boundary of a granule and move together into the granule itself. 
Fig.[5]shows in the top row granules of the visible solar surface in fields of view of 5600 
X 5600 km as observed with the SST and the Sunrise ballon-borne telescope. Events 
of granular lanes are marked with arrows. Simultaneously recorded Doppler maps (see 
ISteiner et al.. 2010) show that material flows in the upward direction within the bright 
granular lane with speeds of up to 1.8 kms~^ Most of the area behind the bright rim, 
including the dark lane, harbors upflows as well but less strong, and sometimes even 
downflows. Close to the border of the granule and within the adjacent intergranular 
lane, the plasma flows in the downward direction with speeds up to 1.2 kms~^ 









IIT73T 







Figure 5. Top row: Granules of the visible solar surface in fields of view of ap- 
proximately 5.6 X 5.6 Mm as observed with the SST (1st panel, courtesy L. Rouppe 
van der Voort) and the Sunrise ballon-borne telescope (2nd and 3rd panel). The 
threadlike, bright/dark stripes, marked with arrows, develop at the edge of a granule 
and move into the granule itself. Bottom row: Equivalent sections from an MHD 
simulation with the C05B0LD code. The arrows mar k strikingly similar t hreadlike 
objects, which prove to be vortex tubes. Adapted from lSteiner et al.l (l2010h . 



The bottom row of Fig. [5] shows equivalent sections from a computer simulation. 
The arrows mark strikingly similar threadlike objects as are observed. Fro r n cros s sec- 
tions through the computational domain of the simulation, ISteiner et all ( 20101) con- 
clude that these granular lanes are the visible signature of horizontally oriented vortex 
tubes. The axis of the vortex tube typically coincides with the trailing dark edge. Di- 
rectly above the vortex tube, the flow assumes transonic speeds, roughly parallel to the 
solar surface. There, gas pressure and temperature are low, which opens a relatively 
transparent view to the cool interior of the vortex tube causing the dark edge. Rather 
than a proper movement, it is the shape and size of the vortex tube which changes in 
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Figure 6. Rise and fall of a vortex arc. The plots display the swirling strength 
(green volume rendering) and the optical surface (yellow) at three different times 
(labels are in mi nutes). The size of the box shown is 1.5 x 1.5 x 0.8 Mm^. Credit: 
iMoU et al] (1201 li) reproduced with permission © ESO. 



time and causes the horizontal displacement of the bright lane and the expansion of the 
trailing darkish area. In the simulation, we see occasioi ial magnetic fi eld intensification 
at locations of vortex tubes — one such case is shown in lSteiner et all (1201 ll) . Whenever 
magnetic fields are present, they tend to point in the horizontal direction near or above 

the Tc - 1 surface. 

lYurchvshvn etal.1 (120111) observe that inter granular jets, origi nating in the inter- 
granular space surrounding individual granules (Good e et al.ll2010l) . tend to be associ- 
ated with the formation and evolution of granular lanes. They speculate that the inter- 
granular jets may result from the interaction of the turbulent small-scale fields associ- 
ated with the vortex tube with preexisting larger-scale fields in the intergranular lanes. 
Since we found in Fig.|4]a singular case where twin whirlpools beneath the surface were 
an integral part of a granular lane, it cannot be excluded that sinkholes are an additional 
important ingredient for the production of the intergranular jets. 



2.5. Small-Scale Vortices in Simulations of Solar Surface Convection 

A synoptical view of sm all-scale vortices in simulations of solar surface convection was 
provided by IMoU et all (201 1 ). These authors define the "swirling strength" based on 
an eigenanalysis of the velocity gradient tensor. It allows them to pick high-vorticity 
regions in which the plasma is swirling. The swirling regions form an unsteady network 
of highly tangled filaments, some of which protrude above the optical surface. 

Near the optical surface, vertically oriented swirls are preferentially located within 
the intergranular lanes, where cooled fluid is sinking down in a turbulent fashion. Hor- 
izontal swirls are predominant at the edges of the granules. Above the optical surface, 
the three-dimensional structure of swirls is manifold, but often they appear in the form 
of bent and arc-shaped filaments as shown in Fig.O A similar type of str ucture was also 
report ed by Stein & Nordlund (1998) and by Muthsam et all d2010.) . Kitiashvili et al. 



(120121) discuss vortex tubes that reach into the chromosphere and interact with magnetic 
fields in their 3D simulations and they conclude that magnetized vortex tubes would ef- 
fectively transfer energy and momentum from the photosphere to the chromosphere. 

From Fig. |6l one can estimate that the arc-shaped swirl has a diameter of about 
30 km, which is close to the grid size of 10 km of the numerical simulation. One could 
therefore speculate that in reality the swirling filaments would be even thinner. It would 
be extremely difficult to detect them with present day solar telescopes. Doppler-shift 
measurements have a depth resolution, limited by the formation width of the spectral 
line or a combination of spectral lines, which is definitively larger than 30 km. 
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3. The Horizontal Magnetic Field of Quiet Sun Regions 

A fundamental achieve ment attained with Hinode's Solar Optical Telescope/Spectro- 
Polarimeter (SOT/SP) dSakuraillloOSh was the discovery that in the mean, the mag- 



neti c field over wide areas of the quiet Sun photosphere is predomin antly horizon- 
tal (iLites et all l2007l. l2008l : lOrozco Suarez etalJ l2007l : iDanilovic et aO i2010). From 



the analysis of time-averaged deep mode Stokes spectra with an effective integration 
time of 67 . 2 s an d a noise level in the polarization continuum of about 2.9 x 10"^ Ic, 
Lites et ahl (l2008h came to the conclusion that the spatially averaged horizontal appar- 



ent flux density, as derived from wavelength-integrated measurements of the Zeeman- 
induced linear polarization, is 55 Mxcm"^, five times more than the corresponding av- 
erage vertical apparent flux density of 11 Mxcm^^. This result, however, remained not 
undisputed. 



3.1. Critiques of the Horizontal Field Measurements 



Martinez Gonzalez et al.l ( 20081) analyzed spectro-polarimetric data of the Fe i lines at 



1.5648 /im and 1.5652 /im, recorded with the Tenerife Infrared Polarimeter (TIP) at the 
Vacuum Tower Telescope (VTT, Observatorio del Teide), of very quiet Sun regions at 
various heliocentric distances on the solar disk from ju = 1.0 to = 0.28. They re- 
strict the analysis to profiles with a degree of polarization, yjQ^ + + larger than 
4 X 10"^ Ic- They found that the circular and linear polarization amplitudes do not have 
any clear dependence on the heliocentric angle. This result points to an isotropic distri- 
bution of magnetic fields and is against any field topology with a preferred orientation 
with in the field-of-view. 

|Beck & Rezaei ('2009') observed a quiet Sun region at disk center with TIP® VTT, 
using the same lines as Martinez Gonzale z et a l. (2008). The noise level of the polariza- 
tion signal was 2 x 10"^ Ic, where they set a threshold of 1 x 10~^ Ic for the polarization 
degree during the analysis. They find that the total magnetic flux contained in the more 
inclined to horizontal fields (y > 45°) is about two times smaller than that of the less 
inclined fields. Ho wever, they do not conclude that this result would necessarily con- 



tradict the results of iLites et all because of the different formation heights of the lines at 
1.56 um and at 630.2 nm. In fact, simulations (Schiissler & Vogler 2008 ; Steiner et all 
20081 : IDanilovic et al.ll2010l) show a strong height dependence of the horizontal mag- 



netic field — while the vertical component may dominate in the deep photosphere, the 
horizontal c omponent becornes imp ortant in the upper photosphere and lower chromo- 
sphere only. iBeck & Rezaeil (l2009l) also emphasize that the pixel to pixel variation of 
the thermal structure of the atmosphere should be properly taken into account during the 
analysis of the polarization signal. This is not taken into account when using a global 
(pixel independent) calibration curve for the derivation of magnetic field strengths as 
was done by ILites et al.1 



lAsensio RamosI (l2009h carried out a Bayes ian analysis of the 'normal mode' data 



(4.8 s exposure time, CTnoise * 1.2 x 10 ^Ic) of iLites et al.l (l2008h . He points out that 



the noise level present in the Hinode SOT/SP observations induces a substantial loss of 
information for constraining the angular distribution of the magnetic field. However, 
the results indicate that the field of pixels with small polarimetric signals has a quasi- 
isotropic distribution. He also finds that the magnetic field strength in the internetwork 
region is clearly in the hectogauss regime with 95% confidence. 
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Table 2. Angular distribution of the internetwork magnetic field as determined by 
various authors. The entries in the column 'angular distribution' must be consid- 
ered an over-all summary, the details being more complex. A pure isotropic angular 
distribution would lead to (Bapp)/(Bapp) = 7r/2. The entries in the lower part of the 
table refer to numerical simulations, for which the ratio (Bapp)/(^app) was derived 
from synthesized Stokes profiles and application of a PSF corresponding to that of 
Hinode/SOT. Numbers in parentheses without PSF. 



no. 


authors 


instrument/ 
simulation 


line 
[nm] 


angular 
distribution 




1 


Lites et al. f2007. 2008) 


SOT/SP 


630 


predominantly horizontal 


5 


2 


Orozco Suarez et al. (2007) 


SOT/SP 


630 


predominantly horizontal 


2.1 


3 


Martinez Gonzalez et ah (2008) 


VTT/TIP 


1560 


isotropic distribution 




4 


Beck & Rezaei (2009) 


VTT/TIP 


1560 


predominantly vertical 


0.42 


5 


Asensio Ramos (2009) 


SOT/SP 


630 


isotropic for weak fields 




6 


Danilovic et al. (2010) 


SOT/SP 


630 


predominantly horizontal 


5.8 


7 


Stenflo (2010) 


SOT/SP 


630 


predominantly vertical 




8 


Ishikawa & Tsuneta (2011) 


SOT/SP 


630 


predominantly vertical 


0.86 


9 


Borrero & Kobel (20JT) 


SOT/SP 


630 


undeterminable 




10 


Borrero & Kobel (2012) 


SOT/SP 


630 


non-isotropic 




11 
12 


Steiner et al. (2008) 
Danilovic et al. (2010) 


h20 
vlO 

Cmf=3 

C+5ver 


630 
630 
630 
630 


predominantly hor- 
izontal 

predominantly hor- 
izontal 


4.3 (2.8) 
1.6(1.5) 
9.8 (3.5) 
4.2 (2.6) 



Y et another way of analyzing the Hinode SOT/SP data was advanced by 'S tenflol 
(l2010h . Ae carried out a thorough analysis of the ratio of the amplitudes of the blue 
wing of Stokes V of Fe i AA 630.15 and 630.25 nm and finds a "magnetic dichotomy" 
with two distinct populations, representing strong (kG) and weak fields. With regard 
to the inclination he finds that the angular distribution is extremely peaked around the 
vertical direction for the largest flux densities, but gradually broadens for smaller flux 
densities, to become asymptotically isotropic at zero flux density. This means that the 
m agnetic field has a predo min antly vertical orientation, q uite opposite to the findings 
of lLites et all(l2007ll2008l) and lOrozco Suarez et all (120071) . We come back to Stenflo's 
met hod in Sect. [Ml 

Ishikawa & Tsunet^ ( 201 lb find, based on Hinode SOT/SP data, a clear positional 



association between the vertical and the horizontal magnetic fields in the internetwork 
region. Essentially, all of the horizontal magnetic patches are associated with vertical 
magnetic patches and half of the vertical magnetic patches are associated with horizon- 
tal magnetic patches. This points to small-scale magnetic loops with bipolar footpoints, 
and to canopy fields of magnetic flux concentrations as an important source of the hor- 
izontal magnetic field. In fact, ti ny emerging loops of g ranul ar and sub-granular scale 
were o bserved with SOT/S P by llshikawa et al.l (l2010l) and iMartmez Gonzalez et al. 



(I2010I) . ISteiner et all (120081) show an example of an emerging sub-granular horizontal 
field patch framed by opposite polarity footpoints of vertical field from their simula- 
tions. They describe it as a consequen ce of convective ove r shooti ng, which transports 
horizontal fields into the photosphere. ISchtissler & Vogle^ ( 2008 ) see the loopy struc- 
ture of the small scale magnetic field in the photosphere as a direct outcome of the 
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surface dynam o which operates in the surface layers of the convection zone (see also 



Steinej (1201 Ol) for a sketch explaining why such a loopy structure leads to an increasing 



dominance o f the horizontal fie ld component with increasing height in the atmosphere). 
Prior to this, iLites et al] d 19961) report isolated, predominantly horizontal magnetic flux 
patches of l"-2" and smaller size, often occurring between regions of weak opposite 
polarity Stokes V profiles, which they suggested to be the emergence of small, concen- 
trate d loops of magnetic flux. 

Ilshikawa & Tsunetal ( 201lh measure for an internetwork region a mean longitu- 
dinal (vertical) magnetic flux density of 8.3 Mxcm~^ and a mean transverse (hori- 
zontal) flux density of 7.1 Mxcm~^, thus a clear dominance of t he vertical fields . 
Note that while the former number is close to the value obtained bv lLites etaLl(l2007L 
2008h andio rozco Suarez et all (120070 . the latter is about 8 times smaller. One rea- 



son for this discrepancy is the stringent noise criteria that was applied to the data by 



Ishikawa & Tsunetal (120111) . Since a transversal field produces a much weaker signal 



in linear polarization than a longitudinal field of equal strength in circular polarization, 
there is a bias towards longitudinal fields when applying a common noise level and 
discarding any signal below it. Correspondingly, Ilshikawa & Tsunetal (1201 lb find only 
10. 1% of a supergranular cell sized subfield having sufficient linear polarization for en- 
tering the analysis, while 43.8% have sufficient circular polarization. On the other hand 
applying no noise level gives preference to transversal fields because pure noise would 
produce spurious large amounts of transversal field. We come back to the problem of 
selection effects in Sect. 13.31 

Table |2]provides an overview of the above discussed various attempts to determine 
the horizontal vs. the vertical mean flux density (or the angular distribution) of the 
internetwork magnetic field. 



3.2. The Photon-Noise Problem 

The main reason for the disparate results highlighted in Table |2] is due to the different 
sensitivity of linear and circular polarization to magnetic fields, in combination with 
the finite sensitivity of the mea surements due to photon noise. The problem was nicely 



demonstrated and discussed by iBorrero & Kobell (1201 ll) . In order to test if it was actu- 
ally possible to accurately retrieve the magnetic field vector employing the Fe i line pair 
at 630 nm in regions wit h very low polarization signals such as in internetwork regions. 



IBorrero & Kobell ( 201 lb carried out several Monte-Carlo simulations with synthetic 



data. They start from a map of magnetic field vectors and associated atmospheric pa- 
rameters, which they obtained in the first place from a Stokes inversion of real SOT/SP 
data. Then, they arbitrarily set the transversal magnetic field to zero, synthesize the 
Stokes profiles from this map and add different amounts of noise. Then they invert 
these data back to magnetic field vectors as if they were actually observational data. 
Clearly, the subsequent inversion should result in zero transversal field. However, this 
is not the case as is shown in Fig. |7] 

The four panels show the results for different amounts of noise that was added to 
the synthesized profiles — from 1.0 x 10"^ 4 (top left) to 1.0 x 10"^ Ic (bottom right). 
2.8 X 10~^ Ir cor responds to the noise level of the least noisy SOT/SP data used by 
Lites et al. I (12001 12OO8). The inversion procedure asks for a threshold of the signal 



to noise ratio. When either Stokes Q, U, or V are above this noise level, the pixel is 
considered for inversion. This ratio was taken 3, 4.5, and 6, but the different curves 
corresponding to these values are almost identical. The curves are histograms of the 
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Figure 7. Histograms of the transverse component of the magnetic field, B^, 
which is purely induced by different levels of noise, cts, indicated in the top left of 
each panel. The different curves correspond to different signal-to-noise threshold for 
a pixel to be considered for inversion. The vertical black arrow indicates the 5-Dirac 
distribution of the true, original transverse component. The vertical, dashed lines 
indicate the center of gravity of the histograms. Credit: iBorrero & Kobell (1201 lb 
reproduced with permission © ESO. 



transversal magnetic flux density and should really be 5-Dirac distributions centered 
at zero if the inversion was correct. However, the center of gravity of the histograms 



(indicated by the vertical dashed lines) is 80-90 G for a noise level of 1.0 x 10 ^ /, 



c 



and close to 55 G for the noise level of 2.8 x 10"^ Ic, which is, coin cidence or not, the 
value for the mean transversal apparent flux density determined by iLites et all (2007, 
20081 ). Even at a hypothetical noise level of 1.0 x 10~^ 7^, the retrieved histogram of the 



transversal field still has a center of gravity of 10 G and peaks at about 6 G instead of 
OG. 

Here is a note in place. The signal to noise limits used here refer to single wave- 
length points, not to integral quantities like 2tot = j^' QW j L j^' ■ Thus, even 
though the probability that photon noise produces a signal above 3o", at a given wave- 
length position, is only 0.3%, it is almost guaranteed that one of the 112 wavelength 
points in either Stokes Q, U, or V surpasses this limit and therefore the full profiles 
enter the analysis. Thus, even with a signal to noise ratio of 5, there are still ample pro- 
files that enter the analysis with linear polarization that is interpreted as real although 
it would be strictly zero without the addition of noise. Also, most of these profiles are 
chosen because at least one of the wavelength points in Stokes-V surpasses the signal- 
to-noise limit and none in Stokes Q or U. This then leads for example in the case of 
cr - 2.8 X 10"^ /c to a residual horizontal field strength of 55 G, much more than the 
corresponding vertical field strength, which in this case is around 20 G and therefore. 



16 



O. Steiner and R. Rezaei 



the field would be interpreted as predominantly horizontal although it is strictly verti- 
cal in reality. This devastating result sheds serious doubts on previous determinations 
of the angular distribution of internetwork magnetic fields based on Stokes inversion 
methods. 

In defense of the work of i Lites et all ( 20071 12008 ) it should be noted however 



that they worked exclusively with integral quantities, Vtot and Ltot, which should be 
much less affected by noise. However this forced them to use a previously determined 
calibration curve for the conversion of total linear and total circular polarization into 
field strength. This calibration was derived from a given atmosphere and therefore it 
cannot take pixel to pixel varia tion of the atmo s pheric paramet ers into acc ount, which 
was found to be important by iBeck & Rezaeil ( 20091) . Also, iLites et all did not use 



Stokes Q and U separately, but tried to determine the "preferred-frame azimuth" for 
each pixel so that Stokes U vanishes, enhancing the signal in Stokes Q. This procedure 
should enhance the signal to noise in l inear polariz ation but is of course itself subject to 



noise. iDanilovic et all (12010 ) applied ILites et aLl 's method to white noise with a level 



of cr = 8 X 10"^ Ic and obtained a mean spurious transversal magnetic field strength of 
36 G, w hich is abo ut half of what is shown for a similar noise level in Fig. |7J indicat- 
ing that'Lites et al.''s method produces considerably less spurious horizontal field than 
straightforward Stokes inversion does. 

3.3. Selection Effects 

There are different ways of selecting pixels as having good enough signal for a trustwor- 
thy determination of the magnetic field vector, all leading to biases in one or the other 
way, again due to the largely different sensitivity of the circular and linear polarization 
to magnetic fields. 

Selecting pixels only which have either Stokes V or Q or U above the noise level 
gives a strong bias towards horizontal fields. This is because most pixels will be selected 
because they have a Stokes V above the noise level (because of the high sensitivity of 
Stokes V to longitudinal fields) but these pixels are dominated by noise in Stokes Q and 
U, which then will yield substantial amounts of false transversal field (because of the 
low sensitivity of linear polarization to transversal fields). 

One could think that the obvious remedy consists in selecting pixels only that 
have signals in Stokes V and Q or U (V A{Qv U)) above the noise level. However, this 
again gives preference to horizontal fields because in this case only pixels with strong 
transversal fields are chosen, while the many pixels with Stokes V above the noise level 
and Q or U below it are discarded. 

One way (which has not been tried out so far) to circumvent these problems would 
consist in setting limits in the real physical space of field strengths, instead of in the 
Stokes space. Thus, we demand that either By or Bj^ be larger or equal to a given 
field-strength limit Bum- Then, 

Viim ^ CcBiiyn and Qiim ^ B^;^ , (2) 

where we now work without loss of generality in the preferred-frame azimuth where 
U vanishes. Cc and c/ are calibration constants that could possibly vary form pixel to 
pixel. Then 

^lim „ /Glim 1 ; ,^^ 
= Biiyn = ^ — ^ = — V^CTnoise , (3) 
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where we now demand that Stokes Q be ncTnoise above the noise level if it is to be used 
for determining the transversal field. This then leads to the selection criterium 

Cc 

Q > MCTnoise or V > — ^|na~ . (4) 

Note that the criterium for Stokes V is much more stringent than that for Q, which is 
expression of setting equal lower field-strength limits to the transversal and longitudinal 
field. Still, we expect many pixels to satisfy the stringent criterium for Stokes V but 
still not the less stringent criterium for Stokes Q. In this case we propose to discard 
Q, assuming the field to be strictly vertical. On the other hand in the rare case where 
Q satisfy the criterium but not V, we consider the field as strictly horizontal. In case 
where both Q and V satisfy their individual criterium, the magnetic field vector can be 
determined. 

These criteria should remedy the one-sided preference for horizontal fields but 
still has its own problems. It can be expected to give some preference to strictly vertical 
and strictly horizontal fields. Also, large amounts of pixels that would have a V signal 
above ncTnoise are discarded for inversion because they do not fulfill the more stringent 
criterium for V, Eq. |4l Therefore, for only a fraction of a given polarization map the 
magnetic field could be retrieved, making it difficult to reliably determine mean values 
of transversal and longitudinal flux over the full map. 

3.4. Stokes-F Amplitude Ratios 

There are, roughly speaking, two ways of retrieving the magnetic field vector from po- 
larimetric data. With the Stokes inversion procedure one starts with a given atmospheric 
model including a magnetic field, computes the Stokes profiles of one or several spec- 
tral lines by solving the Unno-Rachkowsky equations of polarized radiative transfer in 
a forward fashion, and compares them to the observed ones. Then, the magnetic field 
vector, thermal atmospheric structure, line-of-sight velocity, magnetic filling factor, and 
possibly the micro and macro turbulence parameters and stray-light contributions are 
iteratively adjusted until a satisfactory agreement between synthesized and observed 
profiles is achieve d. The first such automatic Stokes inversion code was described in 
iKeller et al.' ('1990'). This procedure is repeated pixel by pixel of a polarimetric map, so 
that the atmospheric parameters are allowed to vary, say from the interior of a granule 
to the intergranular lane. Stokes inversion procedures were chosen by the authors of 
the second, fourth, fifth, and ninth row of Table |2] In a sense, it can be considered 
the classical way of atmospheric modeling as is performed in deriving standard stellar 
model atmospheres. 

Many solar physicist distrust Stokes inversion procedures because of the many free 
parameters that need to be adjusted (which renders the solution non-unique), the poten- 
tial vulnerability to nois e, and the re l atively complex programs, which give the impres- 
sion of black boxes (see iRuiz Cobol (l2007b for replies to such concerns). Instead, they 
take resort in simple procedures derived from first principles of spectro-polarimetric 
line formation theory. This approach may allow them to work with spectrally integrated 
quantities, like Vtot and Ltot mentioned in Sect. 13. 2[ which may be less susceptible to 
noise. However, such basic receipts invariably end up in using a calibration curve for 
converting Stokes signals to magnetic field strengths. Usually, the calibration applies 
globally so that no pixel to pixel variation in the thermodynamic variables is taken into 
account. Such procedures were chosen by the authors of the first and sixth to eighth 
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row of Table [21 The synthetic data of rows 11 and 12 were inverted with lLites et alt s 
method. Of course, also this second, more fundamental way of deriving the magnetic 
field vector can be considered a Stokes inversion procedure — one is tempted to call it a 
"poor man's inversion procedure". Its advantage is that it immediately derives from the 
observed data and its model dependency is less complex. 

One method on this second route for retrieving the true magnetic field strength of 
quiet Sun regions consists in considering the Stokes-V amplitude ratio of two spectral 
lines. Traditionally, this was done with the two lines Fe i AA 525.022 and 524.706 nm, 
which yield the so called magnetic line ratio (Stenflo lQVSi) . The two lines have eff"ective 
Lande factors 3.0 and 2.0, respectively, but are otherwise identical in terms of line 
strength and excitation potential and belong to the same atomic multiplet. Therefore, 
the two lines behave virtually identical with regard to formation height and thermal 
response but they vary in the sensitivity to magnetic fields. 

In the weak field regime in which the Zeeman splitting is much smaller than the 
width of the spectral line (at most a few hundred Gauss), 

V{A) K , (5) 

oA 

where V{A) is the Stokes V profile, the difference between left and right circularly 
polarized intensity and I{A) the Stokes / profile, the regular, natural light intensity (ab- 
sorption profile), g is the effective Lande factor. Equation Q simply derives from the 
leading term of the Taylo r expansion o f the two Zeeman components expressed in terms 
of the / profile (see e.g. IStenflol[l994l . chapter 12.2). Then the Stokes V ratio of two 



lines 1 and 2 is 

V]_ _ gi dh/dA ^ gi_ 

Vi gi dh/dA g2 ' 
which is in particular also valid for the wavelength position where Stokes V assumes its 
peak value. The second relation of Eq. Q is well satisfied only when lines 1 and 2 are 
virtually identical and form under very similar conditions as is the case for the line pair 
Fe I AA 525.022 and 524.706 nm. As long as the magnetic field is intrinsically weak, 
the Stokes-V amplitude-ratio is gi/g2- But with increasing field strength, the Stokes-V 
amplitude starts to saturate, first the line with the larger g-factor while the V amplitude 
of the line with the smaller g-factor continues to grow and Eq. Q breaks down because 
higher order terms in the above mentioned Taylor expansion come into play. At field 
strength in the order of 1 kG, the two amplitudes are of similar size, meaning that the 
amplitude ratio approaches V1/V2 ~ 1. Further increase of the magnetic field strength 
merely increases the Zeeman splitting between the blue and the red lobe of the Stokes 
profile, linearly with field strength. Thus, the line-ratio method provides an easy and 
comprehensible tool for differentiating between weak and strong fields in the regime 
where full Zeeman splitting of the lines has not yet taken place, making it possible to 
dispense with detailed radiation transfer and modeling. It is important to notice that 
the line-ratio method also works when only a small fraction of the spatial resolution 
element is occupied by magnetic fields (small filling factor) so that the method yields 
information about sub-resolution magnetic fields, i.e., the intrinsic field strength. 

Unfortunately, Hinode SOT/SP does not detect the ideal line pair Fe i AA 525.022 
and 524.706 nm but instead the line pair Fe 1 AA 630.151 and 630.250 nm with respec- 
tive Lande factors 1.667 and 2.5, for which the second relation of Eq. dSjl is less well sat- 
isfied because the ratio of Stokes-/ derivatives does not cancel ou t but introduces an in- 
tricate dependence on temperature and velocity in the atmosphere dKhomenko & CoUadosI 
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Figure 8. Scatter plot of the blue lobe Stokes-V amplitude of Fe i 630.250 nm vs. 
the corresponding amplitude of Fe i 630.1 51 nm from the lef tmost third of the 'nor- 
mal mode' map shown in Figs. 1 and 2 of iLites et al.l (l2008h . The dashed line with 
slope i = 1 .66 represents the regression relation that would be expected for intrinsi- 
cally weak fields. Note the two populations of points: (1) The dominant population 
that follows a slope of about 0.65^ for small polarizations. (2) A secondary pop- 
ulation of Stokes Y amplitudes below about 2%, which closely follows the dashed 
line. The me thod for determining the Stokes- V amplitudes was the same as used by 
IStenflol (120101) (Fig. 8) for a different ('deep mode') data set. 



2007h . Therefore it did not seem advisable to apply the line-ratio method to SOT/SP 



data. Nevertheless, Stenflo (2010 ) gave it a try, which led him to a remarkable discovery 
which is displayed in Fig. |8] 

Fig.[8]shows a scatter plot of the blue lobe Stokes-V amplitude of Fe i 630.250 nm 
vs. the corresponding amplitude of Fe i 630.151 nm for a SOT/SP "normal mode" data 
set recor ded at quiet Sun di sk-center, viz., the leftmost 100" of the map shown in Figs. 1 
and 2 in Lites et al. ("2008^. Fig. [8] is from a different dataset than, but very similar to 
Fig. 8 of [Stenflo (2010.) and th e Stok es- V amplitudes were determined with the same 
method as described in Stenflol (l2010h . The dashed line with slope s = 1.66 represents 
the regression relation that would be expected for intrinsically weak fields. It is close 
to the ratio of the respective Lande factors, ^630.25/^630.15 - 2.5/1.667 = 1.5 but not 
exactly because the ratio of the derivatives in Eq. ^ is not unity but rather 1.1 at the 
wavelength position of maximum Stokes V amplitude of the blue lobe. The unexpected 
result is that there exists two populations of points: (1) The dominant population that 
follows a slope that is about 0.65s for small polarizations, decreasing to about 0.55s 
for larger values. (2) A secondary population that is significant only for Stokes V 
amplitudes below about 2%, and which closely follows the dashed hne. Even at the 
1 % polarization level, the two components can be clearly separ ated when plo tting the 
histogram of the points in a bin around this polarization level (IStenfloll2010l Fig. 9). 
Reliable decompositions are found down to V630.15 > 0.6%. 
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Figure 9. Top row: Scatter plot of the Stokes-V amplitude of the blue lobe of Fe i 
630.250 nm vs. the corresponding amplitude of Fe i 630.151 nm of synthesized V 
profiles from the C05B0LD MHD simulation v50. The scatter to the left is from 
pixels at full spatial resolution, the scatter to the right from the same data after appli- 
cation of the SOT PSF The vertical lines indicate the polarization levels which the 
histograms in the bottom row correspond to. Bottom row: Histograms of the points 
that fall within the range 0.8 < Veso.is < 1-0 (black) and 1.8 < Veso.is < 2.0 (gray). 



Fig- m has a simple and elegant explanation in terms of the Stokes-V line-ratio the- 
ory. Clearly, the points of the second population must correspond to intrinsically weak 
fields because their Stokes-V amplitude ratios are in close agreement with the weak field 
limit (dashed regression relation). The points of the first population must correspond 
to kG fields because their Stokes-V amplitude ratios is approaching, and even dropping 
below, one, which points to strong fields as explained above. But how is it possible that 
the strong field component exists not only for large polarization amplitudes, which may 
stem from network elements, but extends down to very small polarization levels? One 
would rather expect the scatter plot to display a S-shape — the inner part being occupied 
by the weak field component with slope s ~ 1.66 and the outer part by the strong field 
component with slope s ~ 1.1. The answer to this question is filling factor. There must 
exist ample amounts of kG-fields of sub-resolution size (small magnetic filling factor), 
which give only weak polarization amplitudes but nevertheless betray themselves as 
being intrinsically strong because of a small Stokes-V amplitude ratio. Thus, pixels 
with polarization levels around and below 1% must harbor kO fiux concentrations of 
very small size and in fact, Stenfio (2011) deduces from these results kG fiux tubes 
with diameters down to 10 km. As elegant and obvious this explanation of Fig.[8]may 
appear, caution is indicated as is demonstrated in the following. 
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3.5. Stokes- y Amplitude Ratios from Simulations 

Fig. |9] shows the inner part of a scatter plot of the same quantities as shown in Fig.[8]but 
of Stokes-V profiles that were synthesized from the C05B0LD MHD-simulaion v50 
referred to in Sect. 12.31 The top left panel shows the scatter plot from the data at full 
spatial resolution, the top right panel after a pphcation of the Hinode/SOT point spread 
function (PSF) of'Wedemeyer-Bohin| (12008') to the full resolution data. Here, we do not 
use Stenflo's fitting method for determining the Stokes- V amplitudes but take the actual 
amplitude, discarding abnormal Stokes-V profiles. The bottom panels show the corre- 
sponding histograms of the points that fall within the range 0.8 < ^630.15 < 1-0 (black) 
and 1.8 < V^^o.is ^ 2.0 (gray), again, once from the full resolution data (left) and 
once from the data after application of the SOT PSF (right). We see immediately that, 
like in the real, observed data, there are two populations of points: (1) The dominant 
population that follows a slope that is about 0.6s. (2) A secondary population that is 
significant only for Stokes V amplitudes below roughly 2%, and which approximately 
follows a regression line with slope s = 1.66. For the full resolution data, even at the 
1% polarization level, the two components can be clearly separated as can be seen from 
the two distinct humps in the corresponding histogram. 
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Figure 10. Left: Scatter plot of the blue lobe Stokes-V amplitude of Fe i 
630.250 nm vs. the corresponding amplitude of Fe i 630.151 nm of synthesized V 
profiles from the C05B0LD MHD simulation h50. The scatter is from pixels at full 
spatial resolution. Right: Corresponding histogram of the points that fall within the 
range 0.8 < Vao.is < 1-0 (black) and 1.8 < Veso.is < 2.0 (gray). 




A scatter plot from a snapshot of the simulation h50 (Fig. [TOb shows a similar 
behavior with the difference that both polarities are equally represented and that the 
weak field population is now the major one. This is because in this simulation the field 
is weaker and more turbulent than in simulation v50 and has no preferred polarity from 
the beginning. 

On the background of the explanation for the observed scatter plot in terms of 
intrinsically weak and strong fields and in terms of the magnetic filling factor, the result 
from the synthesized data is very surprising. Since the simulations cannot harbor any 
sub-resolution magnetic elements because there cannot be any structure smaller than 
the resolution limit given by the computational grid, the magnetic area filling factor for 
a single computational cell is always one, irrespective of the field strength. Thus, the 
argument that pixels of small polarization degree with a Stokes-V line-ratio in the order 
of unity must origin from kG flux concentrations of small filling factor cannot apply to 
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the simulation data. For those, the two populations of points in the scatter plot must 
have a different origin. 



3.6. Interpretation of the Simulation Data 

On the one hand side, it is reassuring that the MHD simulations reproduce the observed 
two populations of points in scatter plots of the Stokes-V amplitudes of the lines Fe i 
AA 630.151 and 630.250 nm. On the other hand, we found in Sect. 13.51 that the two 
populations ca nnot be explai ned in terms of the very attractive and elegant hypothesis 
put forward by IStenflol ( 20101) of a "magnetic dichotomy" with two distinct populations 
representing strong (kG) and weak fields, or "collapsed" and "uncollapsed" fields as he 
expresses it. The question then is, what causes the two populations if not the dichotomy 
of "collapsed" and "uncollapsed" fields? 
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Figure 1 1 . Histograms of the absolute magnetic field strength at three different 
optical depth levels log r = -3 (top), -2 (middle), and -1 (bottom), for the points 
of the full resolution scatter plot of Fig. |9] with a polarization level of 1.0 < Veso.is ^ 
2.0. Dashed curve: Histograms of the main population of points with Stokes-V 
amplitude ration Vg-^Q 25/^530 15 < 1.5. Solid curve: Histogram of the population 
with Stokes-V amplitude ration ^63025/^630 15 - 

Figure [TT] partially answers this question. It shows the histograms of the abso- 
lute magnetic field strength at three different optical depth levels logr = -3 to -1 
for the points of the full resolution scatter-plot of Fig. |9] with a polarization level of 
1.0 < ^630.15 < 2.0. Each panel contains two histograms: one for the main population 
of points with Stokes-V line-ration ^63025/^630 15 - ^-^ (dashed) and one for the pop- 
ulation of points with Stokes-V line-ration ^63025/^63015 - (solid). The second 
population should correspond to intrinsically weak magnetic fields because they cluster 
around a Stokes-V amplitude ratio of 1.6, which is the value in agreement with the weak 
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Figure 12. Top: Location of pixels (white) with a polarization level 1.0 < 
^630.15 ^ 2.0 belonging to the population of pixels with ^63025/^63015 - ^-^ 
(first, main population). Bottom: Location of pixels with a polarization level 
1.0 < Vpo.15 ^ 2.0 belonging to the population with ^63025(^63015 - 
ond, intrinsic weak field population). Background: Continuum intensity at 630 nm. 
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field approximation. In fact, we can see from the liistograms tliat this is true: the field 
strength for this population stays essentially below 200 G with a peak at 10 G. But we 
also see that the first population that should correspond to "collapsed" kG flux concen- 
trations is weak field as well with field strengths mostly below 300 G. There is however 
a distinct difference with respect to the second population: the peak and bulk of the 
histogram shifts to higher field strengths with height in the atmosphere. This means 
that for this population of points, the magnetic field strength increases with height in 
the atmosphere, which is opposite to the usual case. 

When plotting the histogram of the magnetic field inclination with respect to the 
vertical direction (not shown here) we see that the first population shows mainly hor- 
izontal fields in the deep layers, which become inclined by ~ 45° at logr = -2 to 
-3 when the field strength increases. This behavior points toward s cano py magnetic 
fields as the origin of this, formerly called, the "collapsed" (StenfldHoiS) component. 
Canopy fields may stem from magnetic flux concentrations in the photosphere, which 
expand and fan out in the lateral direction with increasing height in the atmosphere 
(ISteinej |2000). The magnetic field along a line of sight that traverses the canopy field 
is weak in the optically deep layers below the canopy but abruptly increases in strength 
and assumes a large inclination angle at the transition from below to within the canopy. 
Another indication that points to canopy fields comes from the histograms of the ver- 
tical velocity, which show mainly upflows at all optical depth levels for the second, 
intrinsically weak field population, but downflows in the deepest layer of the first popu- 
lation. This can be expected because canopy fields are most likely found in the vicinity 
of (vertical) magnetic flux concentrations which in turn are found within intergranular 
lanes, which harbor downflows. Therefore, we expect canopy fields and hence, pixels 
belonging to the first population, to exist at the border between intergranular lanes and 
granules with a tendency of having downflows in the deep layers. 

Figure [12] confirms that this is indeed the case. The top panel shows on the back- 
ground of the continuum intensity at 630 nm the location of pixels that have a polariza- 
tion level 1.0 < ^630.15 ^ 2.0 and belong to the first population with ^53025/^630 15 - 
1.5. Clearly, these pixels are found along the boundaries of granules, where canopy 
fields exist. The bottom panel shows the pixels belonging to the second population with 
^630 25^^630 15 - These pixels occur in patches within granules and are mainly as- 
sociated with horizontal field, which is probably transported into the photosphere by 
instances of granular overshooting. 

From all this we conclude that the two population of pixels with low polarization 
amplitudes found in data of synthesized Stokes- V profiles of magneto-convection sim- 
ulations are not due to a dichotomy between "collapsed" and "uncoUapsed" fields, but 
to a distinct difference between profiles from granule interiors and profiles from the 
boundary region between granules and the intergranular space. Truly strong fields are 
found within the intergranular lane proper but they have a larger polarization degree 
and therefore pertain to the outskirts of the line-ratio scatter-plots. It is very likely, that 
this explanation of the two populations also holds for the observed Stokes-V profiles. 
This then means that there is no need for the i ntroduction of collapse sub-resolution 
flux-concentrations, and that the conclusions of IStenfld (I2OIOD that the inter-network 
magnetic field was predominantly vertical is premature because of its shaky premises. 
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3.7. The Horizontal Field from Numerical Simulations 

From the critiques reviewed in Sects. 13.11 and 13.21 it is clear that the angular distribu- 
tion of the magnetic field in network-cell interiors has not yet been reliably measured. 
It is not even clear if the re exists indeed a predo minance of the horizontal field com- 
ponent. The analysis of iLites et aLl (I2007L |2008|) . which uses and derives least noise 



affected data, points in this direction. On the other hand, numerical MHD simula- 
tions of the solar surface layers have resulted in a predominance of the horizontal fields 
dOrossma nn-Doerth et all 1 19981 : lAbbettll2007l : ISchussler & Vogled l2008l IStemer et al. 
I2OO81 : iDanilovic et al.ll201ol7 particularly in the higher layers of the photosphere and 



low chromosphere. Here, we corroborate these findings with the simulation runs men- 
tioned in Sect. 
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Figure 13. Left: Horizontally averaged strength of various magnetic field com- 
ponents as a function of geometrical height in the model atmosphere of simulation 
h50. Dashed: Vertical component, (|Bvl); Dot-dashed and dotted: Horizontal com- 
ponents <|B.v|> and (IB,!); Solid: Horizontal component Bh - {(BI + By)'''^). Right: 
Histogram for the cosine of the inclination of the magnetic field with respect to the 
vertical direction for different optical depth levels from log r = (broadest distribu- 
tion), to log T - -3 (narrowest distribution). An isotropic distribution would give a 
straight horizontal line. 



Figure [13] shows in the left panel the horizontally averaged strength of various 
magnetic field component as a function of geometrical height in the model atmosphere 
of simulation h50. The dashed curve refers to the absolute value of the vertical compo- 
nent. By, the dotted and dash-dotted curves to the absolute value of the two horizontal 
components Bx and By which stand orthogonal to each other, and the solid curve to the 
horizontal component B^ - (B^ + ^y)''^^- The horizontal components surpass the verti- 
cal one throughout the photosphere from z = to 500 km. They are particularly domi- 
nant at and closely above 500 km height. At z = 500 km <Bh>/<|Bvl>(500km) = 2.9 and 

<S;,,3,>/<|Bverl>(500km) - 1.9. 

The panel on the right hand side of Fig. [13] shows the histogram for the cosine of 
the inclination of the magnetic field with respect to the vertical direction for different 
optical depth levels from logr = (broadest distribution), to logr = -3 (narrowest 
distribution). An isotropic distribution would be represented by a straight horizontal 
line. Instead, the distributions aie peaked at the 90° angle and become narrower with 
increasing height in the atmosphere, because of the increasing dominance of the hori- 
zontal component with height. 
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Figure 14. Horizontally averaged strength of various magnetic field components 
as a function of optical depth of the model atmosphere of simulation h50 (left) and 
the surface-dynamo simulation of lSchiissler & Vogled (12008 ) (right, reproduced with 
permission © ESO). Dashed: Vertical component, {|Bvl); Dash-dotted and dotted: 
Horizontal components and (|Bvl); Solid: Root-mean-square of the horizontal 
component Bh.ims = ({B] + 



The situation is different for simulation v50. Even though, the horizontal compo- 
nents show distinct local maxima close to z = 500 km, the vertical component domi- 
nates for all heights. This is because of the initial condition of a homogeneous vertical 
magnetic field Bo, which conserves the mean vertical field strength above the convec- 
tion zone to Bo for all times. Therefore, the ratio (Bh)/(|fivl) is basically predetermined 
by the initial condition in this case, which may be more representative of a network- 
field patch or a weak plage region. 

Figure [14] shows a comparison between the simulations of I Schussler & Vogler 

( 2008h (right panel) and the present result from simulation h50 (left panel). Again, the 
mean horizontal and vertical components are plotted, this time as a function of con- 
tinuum optical depth and on a logarithmic scale. Despite the very different initial and 
boundary conditions and despite the different codes with which theses simulations were 
carried out, the general trend is similar. In both simulations, the horizontal field com- 
ponent dominates throughout the photosphere. The absolute values of the field strength 
are different though: in case of the simulation h50, it depends on the strength of the 
horizontal field that is advected across the l ower boundary (which is a boundary condi- 
tion), in case of ISchussler & Vogled (l2008l) it depends on the efficiency of the surface 
dynamo, which in turn depends on the Reynolds numbers achieved by the simulation. 



4. Summary and Conclusions 



Vortical flows of mainly vertically directed v orticity hav e been obser ved to exist in the 
deep photosphere (jBrandt et al. 1988; Wans e t al.|ll995l: iB onet et a Lll20()8l: lAttie et al. 
2009| : iBonet et all lgOlO: Balmaceda et al. 2010; Varaas Dominguez et al.ll201 l') and in 



the chromosphere (Wed emeyer-Bohm & Rouppe van der Voort i 200C " Numerical sim- 
ulations show particularly strong vortical flows immediately below the visible optical 
surface Tc = I, beneath (the vertices of) inter-granular lanes. The centripetal force as- 
sociated with these vortical flows is so strong that the gas pressure gradient opposing it 
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leads to a subs tantial reduction of the density within them, which may h ave observable 
consequences (lNordlundlll985l :lKitiashviH et alJ l20TTl : lMoll etal]|201lh . We find that 
there exists no one-to-one connection between the density deficient vortices in the sur- 
face layers of the convection zone and the observable vortical flows at Tc = 1, but one 
can assume the latter to extend and intensify with depth, simply because they are found 
within downdrafts and because of angular momentum conservation in combination with 
the gravitational stratification. This indicates that vortex s tretching is pr o bably an im- 
portant source of vorticity in the surface layers although 'Shelyag et al.| (12011 ) found 
the baroclinic vorticity generation to dominate in the convection zone and at the visible 
surface. 

As was anticipated by iNordlundl (11985 ) and iBrandt et al. disss'), these vortical 
flows can be an important source of hydromagnetic disturbances. Magnetic field that 
is advected into vortices start to rotate. In the photosphere, where the thermal pressure 
generally dominates the magnetic energy density (yS >> 1), this can be considered a 
kinematic process, but vertically extending magnetic fields quickly start to dominate 
with height in the atmosphere. When /3 « I, the rotating magnetic field now acts on 
the plasma via Lorentz force making it to rotate. Therefore, in the photosphere and in 
the chromosphere, the surface of yS = 1 separates two regimes of vorticity generation 
(see also Fig. O. In the regime, where j3 « I vortical flows of vertically directed 
vorticity is generated by vertically extending magnetic flux concentrations which root 
within vortical flows of the deep photosphere and the surface layers of the convection 
zone. Thus, while vorticity is generated by hydrodynamic processes (vortex stretching 
and baroclinic vortex generation) in the deep photosphere and the convection zone, it 
is mediated to the upper photosphere and chromosphere via magnetic fields. One can 
imagine that the formation of rotating ma gnetic flux concentrations may be the source 
of torsional Alfven waves dJess et alJi2009l) . which could potentially create a substan- 
tial Poynting flux in the outward direction, leading to coronal disturbances (but see 
(^outh et al. 2010) for the cutoff frequency of torsional tube waves) and having con- 
sequences for solar- wind accelera tion and coronal heating in regions of open magnetic 
field configurations ( ZirkedfTQQSt) . 

Granular lanes are composed of a leading bright rim and a trailing dark edge, 
which form at the boundary of a granule and move together into the granule itself. 
Virtually every granule harbors one or several granular lanes during its lifetime and 
they are well visible in high resolution broad band filtergrams of the solar surface. They 
are the visib le manifestation of horizontally extending vortex tube found in numerical 
simulations (ISteiner et al.ll20IC)l) . 

A more comprehensive view of vortices in the solar atmosphere reveals that there 
exists not only vertical and horizontal vortices (which are probably easiest to detect) but 
an unsteady network of highly tangled filaments of swirls, some of which protruding 
above the optical surface, for example in the form of arc shaped filaments (iMoU et al.l 
20Tll) . 



Attempts to measure the angula r distribution of the magnetic field in quiet regions 
and internetwork regions of the Sun ('Lites et al."2007',l200Sl; IOrozco Suarez et al. 1I2007L 



Martinez Gonzalez et al. 200S; Beck & Rezaei 2009; Asensio Ramos 2009; Danilovic et al. 
2010l : IStenfloll2010l : IIshikawa & Tsunetall201ll : lBorrero & Kobelll2012D have led to con- 



troversial results. Basic difficulties arise from the fact that in such regions, (i) the 
Zeeman splitting of spectral lines is much smaller than the line width because the 
field strength is mostly weak or spatially not resolved (small filling factor) so that the 
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field strength cannot be determined by simply measuring the Zeeman splitting. In- 
stead, detailed radiation transfer and atmospheric modeling is required for correctly 
interprete polarimetric measurements. Even more problematic for the determination 
of the inclination angle is (ii) that the linear polarization depends in the weak field 
Umit quadratically on the field strength while the circular polarization has a linear 
dependency. Since both polarization states are subject to a common noise level, the 
transversal field component is more severely affected by noise than the longitudinal 
component. Monte Carlo simulations with different noise levels added to synthesized 
Stokes profiles (B orrero & Kobel 2011.) show, that a Stokes-inversion procedure may 
interpret pure noise as transversal fields with a strength of, e.g., 80 G for a noise level 
of cr = 1.0 X 10"^^ Ic- We propose a less biased determination of the mean horizontal to 
vertical field strength by applying a common noise criterion to the field strengths in the 
real physical domain instead of to signals in the Stokes space. 

Establishing a bivariate scatter plot of the blue amplitudes of Stokes V of Fe i AA 
630.15 nr n and 630.25 n m for each pixel of a quiet Sun region, reveals two populations 
of points ( Stenfloll2010h : (i) a population with a slope that is compatible with the weak 
field limit and (ii) a population that seems to be only compatible with kG fields. We 
have analyzed synthesized Stokes- V profiles from an MHD-simulation with the result 
that they too show the same two populations in corresponding scatter plots but can- 
not be explained in terms of intrinsically strong and weak fields. Rather, we find both 
populations stemming from intrinsically weak fields but from different locations. The 
population compatible with the weak field limit typically stems form clusters of pixels 
located within granules. The population that was formerly interpreted as being due to 
intrinsically strong field of small filling factors stems from pixels located at the bound- 
aries between granules and the inter-granular space which harbor (weak) canopy fields 
with a field strength that abruptly increase with height in the atmosphere and show a 
downdraft in the deep photospheric layers (below the canopy field). So far, we have 
not analyzed the details of the radiation transfer that creates under these circumstances 
Stokes-V profiles that behave, regarding their amplitude ratio, like stemming from kG- 
fields but in reality originate from weak fields. The two lines Fe i AA 630.15 nm and 
630.25 nm do not form under the same conditions, in particular not at the same height 
in the atmosphere. Therefore, it can be either that the two lines sample different field 
strengths because of the canopy character of the field, which implies steep field gradi- 
ents, or that they are differently affected by the downflows in the deep layers, creating 
different Stokes-V asymmetries, or that thermal effects play a role. From the present 
analysis we conclude that the two populations are not due to a dichotomy between 
"collapsed" and "uncollapsed" fields but due to a distinct difference in the formation of 
profiles and in the physical conditions at the location of the formation, i.e., the granule 
interiors and at the boundaries of granules, where canopy fields prevail. It is very likely, 
that this explanation also holds for the two populations of actually observed Stokes-V 
profiles. 

Finally, we show, that a new simulation in which convective updrafts across the 
lower boundary advect horizontal fields of 50 G strength into the computational do- 
main, yield a dominance of the horizontal vs. the vertical field throughout the pho- 
tosphere, increasing with height. This is in agr e ement with and cor robor ates previous 
simulation results by lSchiissler & Vogled ( 2008h . lSteiner et al.l ( 2008l) . and lDanilovic et al" 
(12010) . 
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